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ABSTRACT 


Early data from JWST have revealed a bevy of high-redshift galaxy candidates with unexpectedly high stellar masses. An 
immediate concern is the consistency of these candidates with galaxy formation in the standard cosmological model. In the 
ACDM paradigm, the stellar mass (M,) of a galaxy is limited by the available baryonic reservoir of its host dark matter halo. The 
mass function of dark matter halos therefore imposes an absolute upper limit on the number density n(> Mx, z) and stellar mass 
density p,(> Mx, z) of galaxies more massive than M, at any epoch z. Here I show that the most massive galaxy candidates in 
JWST observations at z ~ 7 — 10 lie at the very edge of these limits, indicating an important unresolved issue with the properties 
of galaxies derived from the observations, how galaxies form at early times in ACDM, or within this standard cosmology itself. 
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1 INTRODUCTION 


ACDM-like cosmological models share a similar basic assumption: 
baryons and dark matter are well-mixed at very early times, and 
as baryons collapse into dark matter halos, the maximum amount of 
baryonic material within a halo will be equal to My = fp Mhalo, where 
fo = Qp/Qm is the cosmic baryon fraction. This, in turn bounds the 
total stellar content of a dark matter halo: My (Mhalo) < Mp(Mhalo)- 
I show how this simple relation can be used as a stringent test of 
either cosmological models with minimal assumptions about galaxy 
formation or the reliability of photometric selection and physical 
characterization of high-redshift galaxy candidates. My analysis is 
in many ways similar to Behroozi & Silk (2018), who connected 
cumulative number densities of dark matter halos to high-redshift 
galaxy stellar mass functions (see also Steinhardt et al. 2016), though 
I also consider the maximal cumulative stellar mass density allowed 
in ACDM. The question of the consistency of stellar mass functions 
and the underlying cosmological dark matter halo mass functions has 
become considerably more urgent with the release of the first data 
from JWST, and with it, a swarm of high-redshift galaxy candidates 
(Adams et al. 2023; Atek et al. 2023; Castellano et al. 2022; Donnan 
et al. 2023; Finkelstein et al. 2022; Labbé et al. 2022; Morishita & 
Stiavelli 2022; Naidu et al. 2022; Yan et al. 2023). 


2 ASSUMPTIONS 


I adopt the base ACDM model of Planck Collaboration et al. (2020), 
which assumes no spatial curvature and initial conditions that are 
Gaussian and adiabatic, as the standard cosmological model. I use 
best fit values for cosmological parameters based on the Plik 
TT,TE,EE+lowE+lensing likelihood applied to the full-mission data. 
The relevant parameters and values for this work are the present-day 
Hubble constant, Hy = 67.32 km s7! Mpc7!; the z = 0 density pa- 
rameter for matter, Qm = 0.3158 (which includes baryons, dark 
matter, and non-relativistic neutrinos); the slope of the primordial 
power spectrum of density fluctuations, ns = 0.96605; the rms am- 


© 2023 The Authors 


plitude of the linear matter power spectrum at z = 0 as measured in 
spheres of radius 8 h~! Mpc, og = 0.8120; and the cosmic baryon 
fraction, fe = Qp/Qm = 0.156 (Planck Collaboration et al. 2020). 

With these values, the linear matter power spectrum is specified at 
all times relevant for structure formation. The non-linear density field, 
home to the dark matter halos that host galaxies, must be computed 
numerically. However, a long line of research starting with Press & 
Schechter (1974) has been devoted to connecting the abundance of 
dark matter halos as a function of redshift and mass to the underlying 
linear matter power spectrum. In what follows, I use the Sheth & 
Tormen (1999) dark matter halo mass function dn(M, z)/dM — 
the number of dark matter halos of mass M per unit mass per unit 
comoving volume at redshift z — to compute the comoving number 
density of halos above a given halo mass threshold, 


i dn(M, z) 


n(> Mand = f dM 


(1) 
Mhato dM 


and the comoving mass density in halos more massive than Mpaio, 


Sa dn(M,z 
pm(> Mawd = f amm FED 


(2) 
Mhalo dM 


These translate directly to upper limits on the statistics of galaxies 
through the straightforward assumption that the largest stellar content 
a halo can have given its cosmic allotment of baryons is My,max = 
fo Mhalo. More generally, we may write Mẹ = € fb Mhalo, with e < 1 
being the efficiency of converting baryons into stars. 

The cumulative comoving number density of dark matter halos 
more massive than Mhalo thus sets an upper limit on the comoving 
number density of galaxies more massive than Mx, 


Ngal (> Mx) < Mhaio(> Mx/ fo) - (3) 


Similarly, the cumulative comoving density of collapsed mass sets 
an upper limit on the density of collapsed baryons, pp(> Mhalo) = 
fo Pm(> Mhalo), which in turn strictly bounds the comoving mass 
density of stars contained in halos more massive than Mpalo; 


Px(> Mpaio) < fo Pm(> Malo) » (4) 
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Figure 1. Limits on the abundance of galaxies as a function of redshift. Curves show the relationship between M, and z at fixed cumulative halo abundance 
(left) and fixed pp(> Mhaio), or equivalently fixed peak height v (right). The most extreme L22 galaxy candidates are shown as blue stars, with uncertainties 
indicating 68% intervals (symmetric about the median) of the posterior probability distribution. The existence of a galaxy with Mẹ at redshift z requires that 


such galaxies have a cumulative comoving number density that is at most the number density shown in the left panel, as those galaxies must reside in host halo of 
mass Mhalo = Mx/( fp €). The cumulative comoving number density corresponding to an observed M, will likely be (much) smaller than is indicated here, as 


the curves are placed on the plot by assuming the physically maximal € = 1. For smaller values of €, the curves in each panel move down relative to the points 


by a factor of e (as indicated by black downward-facing arrows). The right panel demonstrates that even for the most conservative assumption of € = 1, the data 
points correspond to very rare peaks in the density field, implying a limited baryonic reservoir that is in tension with the measured stellar masses of the galaxies. 


and the density of stars contained in galaxies above a given Mx, 


Px(> Mx) < fo Pm(> Mx/ fo) - (5) 


3 RESULTS 


The left panel of Figure 1 shows the relationship between the max- 
imal inferred stellar mass for a given Maio, Mx = fo Mhalo (ie. 
assuming the maximal € = 1), and redshift z for fixed cumulative 
comoving halo number densities ranging from 10710 Mpc™3 (light 
gray) to 107? Mpc™3 (yellow). The curves evolve rapidly with red- 
shift, with the maximal stellar mass corresponding to a fixed cu- 
mulative comoving halo number density increasing by two orders 
of magnitude from z = 20 to z = 8. This rapid rise indicates that 
the mass reservoir available for the most massive galaxies increases 
quickly with redshift at fixed halo number density. The two most 
massive high-redshift galaxy candidates from the Labbé et al. (2022, 
hereafter L22) sample, at z ~ 7.5 (My ~ 10!! Mo) and z ~ 9.1 
(My. = 1010.5 Mo), are shown as blue stars. These objects are unex- 
pectedly massive, with stellar content reflective of halos that have cu- 
mulative comoving number densities no higher than ~ 1075-2 Mpc? 
(if € = 1); for € = 0.32 (0.1), the implied number density is ~ 1077 
(io?) Mpc™?. By comparison, the candidates were found in a sur- 
vey of 38 arcmin’, a volume of V ~ 10° Mpc? at each of the redshift 
bins — 7 < z < 8.5 and 8.5 < z < 10 — considered by L22. 

The right panel of Figure 1 recasts the issue in terms of the 
scarcity of systems as measured by cumulative mass density. In 
extended Press-Schechter models, the peak height v(Mbalo.Z) = 
6¢/7(Mhalo, Z) of an object — where ôe ~ 1.7 is the collapse thresh- 
old and o2(Mpalo, zZ) is the variance of the linear density field at 
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redshift z smoothed on a scale containing an average mass of Mhalo 
— is a measure of the fraction of mass in the Universe contained 
in virialized objects more massive than Mpalo at redshift z. Typical 
halos at z have v = 1, which corresponds to 24% of the mass in the 
Universe residing in halos at least that massive; larger values of v 
indicate increasingly massive and therefore rare peaks in the den- 
sity field at that epoch. The comoving baryon density for each peak 
height in the figure is given in the legend; multiplying this number by 
the volume of a survey gives the total amount of baryons contained 
above the mass corresponding to that peak height and redshift. The 
L22 galaxies have peak heights of at least v = 4.5 (assuming € = 1), 
meaning that at most a fraction 6.2 x 1075 of the baryons in the Uni- 
verse are contained in halos massive enough to host these galaxies. 
For reference, v = 4.5 at z = 0 corresponds to Mhalo ¥ 5 X 10!5 Mo. 
Adopting more reasonable efficiencies of € = 0.32 or 0.1 results in 
rarer peaks with v ~ 5.4 or 6.4. 


Figure 2 shows the cumulative stellar mass density reported by 
L22 at z ~ 9 (left) and z ~ 7.5. The data, which come from in- 
dividual massive objects, lie at the extreme of ACDM expectations 
even in the most optimistic scenario: at both redshifts, the measure- 
ments lie at the theoretical limit of p4(> Mx) = fo pm(> Mx/fo), 
implying physically implausible values of e(z ~ 9) = 0.99 and 
e(z = 7.5) = 0.84. When considering the 1 o error (which incorpo- 
rates uncertainties in the stellar mass estimation, Poisson fluctuations, 
and sample variance), the data become marginally consistent with 
the available baryon reservoirs for an efficiency of e(z ~ 9) > 0.57, 
which is likely an unrealistically high value. Assuming a more plausi- 
ble value of e = 0.1 or 0.32 yields a strong discrepancy with ACDM 
expectations at both redshifts even when considering observational 
uncertainties. 
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Figure 2. Stellar mass density limits. The comoving stellar mass density contained within galaxies more massive than M, at z ~ 9.1 (left) and z = 7.5 (right) 
for three values of the assumed conversion efficiency € of a halo’s cosmic allotment of baryons into stars. Only if all available baryons in all halos with enough 
baryons to form the galaxies reported by L22 have indeed been converted into stars by that point — an unrealistic limit — is it possible produce the stellar mass 
density in the highest M, bin at z ~ 9 measured by L22 in a typical volume of a ACDM Universe with the Planck 2020 cosmology. Results are similar at 
z = 7.5. For more realistic values of €, the required baryon reservoir is substantially larger than the theoretical maximum in this cosmology. When considering 
shot noise and sample variance errors (which comprise the plotted uncertainties on the L22 data points in each panel), the measurements are consistent with the 
base ACDM model if € > 0.57, which would still imply incredibly efficient star formation in the high-redshift Universe. 


4 DISCUSSION 


The first glimpse of high-redshift galaxy formation with JWST has re- 
vealed surprisingly massive galaxy candidates at early cosmic times. 
These systems provide a way to test a bedrock property of the ACDM 
model (or, e.g., assumptions in derivations of stellar masses or the 
viability of high-redshift galaxy candidates): the stellar content of ha- 
los should not exceed the available baryonic material in those halos. 
This requirement does not rely on assumptions such as abundance 
matching but rather is simply a statement about the distribution of 
virialized mass in the Universe as a function of redshift and the bary- 
onic reservoirs associated with those virialized halos: galaxies of 
mass M, can only form if halos of mass M,/(€ fọ) have formed. It 
is also more stringent than the requirement that the observed galaxy 
UV luminosity function not exceed the theoretical maximum com- 
ing from a nearly instantaneous (10 Myr) conversion of a halo’s full 
baryonic reservoir into stars (Mason et al. 2022), as it is an integral 
constraint as opposed to a differential one. The massive, high-redshift 
galaxy candidates cataloged in L22 lie at or just beyond the stellar 
mass density constraint in ACDM. 

There are several sources of observational uncertainty that enter 
these results. The flux calibration of NIRCam is continually being 
updated; L22 use calibrations that take into account updated detector 
offsets that are not yet part of the official JWST reduction pipeline 
(see, e.g., Boyer et al. 2022 for examples of this effect and Nardiello 
et al. 2022 for related discussions of empirical point spread function 
modeling for JWST). With NIRCam photometry, a Balmer or 4000 A 
break at z ~ 5 can be mistaken for a Lyman-a break at z = 12 (Zavala 
et al. 2023); the L22 sample was selected to contain both Lyman and 
Balmer breaks, however, and is at low enough redshift (relative to 
z ~ 15 sources) that NIRCam filters can typically exclude z ~ 5 pho- 
tometric solutions. The resulting photometric redshift estimates have 


single, narrow (o;z ~ 0.25) peaks. The masses of the galaxies are 
computed using the median of four methods for fitting the photome- 
try (see L22 for details) and assume a Salpeter (1955) IMF. Different 
assumptions about the photometry (in particular, properties of neb- 
ular emission lines) or IMF could affect the derived stellar masses, 
with the latter being a particularly intriguing possibility. The mass 
of the candidate at z ~ 7.5 was also corrected for the possibility of 
amplification by mild gravitational lensing; this effect is estimated by 
L22 to be 0.15 dex, and the reported mass (and stellar mass density) 
of this object are therefore reduced by this amount to compensate. 
The error bars in Figure 2 include errors in the volume estimates 
coming from both sample variance and Poisson noise, with the lat- 
ter always being dominant in the regime considered here (Trenti & 
Stiavelli 2008; Behroozi & Silk 2018). 


The discrepancy between the observed high-redshift galaxy can- 
didates and ACDM expectations is robust to uncertainties in cos- 
mological parameters in the base ACDM model: the precision on 
each of the relevant parameters is at the < 1% level (Planck Col- 
laboration et al. 2020). Intriguingly, extensions to the base ACDM 
with enhanced values of og and the physical matter density Qh? 
— such as some Early Dark Energy (EDE) models whose aim is to 
resolve the Hubble Tension — predict earlier structure formation and 
a higher abundance of halos at fixed mass at high redshift (Klypin 
et al. 2021), which would enhance the baryonic reservoirs available 
for forming early massive galaxies. Taking the best-fit EDE param- 
eters from Smith et al. (2022), the cumulative comoving baryonic 
density contained in halos more massive than Mhalo = Mx/ fi, for 
the most massive L22 galaxy candidate at z ~ 9.1 is a factor of 
3.3 larger in EDE than in base ACDM, which is non-negligible; the 
L22 data points would then lie at € = 0.72 instead of € = 0.99. 
However, this EDE cosmology is in stronger tension with values of 
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Sg = 0g YQm/0.3 measured at low redshift and predict that the Uni- 
verse is ~ 13 billion years old (as opposed to 13.8 billion years in the 
base ACDM model), which is in moderate tension with the measured 
ages of ultra-faint galaxies and globular clusters (Boylan-Kolchin & 
Weisz 2021). 

At the redshifts studied here, z ~ 7 — 10, the Sheth-Tormen mass 
function overestimates the abundance of massive halos by 20 — 50% 
relative to numerical simulations (Reed et al. 2003; Despali et al. 
2016; Shirasaki et al. 2021; Wang et al. 2022), meaning their true 
abundance at high redshift is likely lower than the Sheth-Tormen pre- 
diction and the constraints derived here are conservative. However, 
the lack of detailed comparisons between theory and simulations at 
high redshifts and high masses points to the importance of continued 
theoretical work in understanding the universality and applicability 
of halo mass function parameterizations in regimes relevant for JWST 
observations (and other forthcoming observatories). 

The tension discussed in this paper is straightforward: the masses 
measured by L22 are only consistent with expectations from the 
standard cosmological model at the reported redshifts if star forma- 
tion in the earliest phases of galaxy formation is incredibly efficient 
(e > 0.57). In the low-redshift Universe, such efficiencies are never 
seen, with € < 0.2 for all galaxies. The theoretical expectation is 
that efficiencies do indeed increase at high redshift (Tacchella et al. 
2018), though € = 0.57 is still highly extreme and likely implausibly 
high. If the explanation of the L22 galaxies is indeed a very high 
star formation efficiency, it implies that the star formation histories 
of such systems must rise steeply with time, following the behavior 
of the baryon reservoirs inside of virialized structures in ACDM. 
The results presented here could also be explained if the stellar IMF 
differs substantially from the assumed Salpeter form, the flux calibra- 
tion of NIRCam changes from the latest post-flight determinations, 
or the volumes currently surveyed turn out to be highly atypical. 

If none of these explanations holds up and these massive galaxies 
are spectroscopically confirmed, they will pose a serious challenge 
for ACDM structure formation with parameters given by Planck Col- 
laboration et al. (2020) because they signify the existence of a larger 
reservoir of collapsed baryons than is possible in this model. Forth- 
coming wider-field JWST surveys, along with JWST spectroscopy 
of massive galaxy candidates, should be able to quickly confirm or 
refute the existence of this tension. Furthermore, the compatibility of 
any additional high-redshift galaxies or galaxy candidates discovered 
in JWST observations with ACDM expectations can be assessed in 
a straightforward way via Fig. 1. If analysis of JWST data continues 
to reveal the presence of strikingly massive galaxies at very early 
cosmic epochs, more exciting surprises lie ahead for the fields of 
galaxy formation and cosmology. 


DATA AVAILABILITY 


Data from L22, including M» estimates and photometric redshifts, 
are available at https://github.com/ivolabbe/red- massive-candidates; 
this paper uses data from sample_revision3_2207.12446.ecsv, 
commit 59fbbfa (from 2023.01.02). All calculations that go into 
the figures in this paper will be made publicly available at https: 
//github.com/mrbk/JWST_MstarDensity. 
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